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PACS 45 . 70 . -n - Granular systems 

PACS 45 . 70 . Cc - Static sandpiles; granular compaction 

PACS 61 . 43 . - j - Disordered solids 

Abstract. - We measure the two-point correlation of free Voronoi volumes in binary disc packings, 
where the packing fraction (?iavg ranges from 0.8175 to 0.8380. We observe short-ranged correlations 
over the whole range of 0avg and anti-correlations for 0avg > 0.8277. The spatial extent of the anti- 
correlation increases with cj>avg while the position of the maximum of the anti-correlation and the 



extent of the positive correlation shrink with (f>.^ 
corresponds to dilatancy onset in this system. 



,g. We conjecture that the onset of anti-correlation 
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' Introduction. — Dry granular matter interacts only 
via clastic and frictional forces, which require particles to 
be in contact; spatially extended interactions like Van der 
Waals forces typically play no role. However, granular 
particles form networks of force chains [T], which implies 

'he existence of local correlations. Lechenault et al. have 
shown that the logarithm of the distribution of free vol- 
umcfQ in granular matter scales in a non-extensive way 
with the cluster size, which implies the existence of corre- 
lations between Voronoi cells [2] . A similar scaling was 
also observed in monodisperse sphere packings l^. These 
experimental observations raise the question how the spa- 
tial extent of such correlations changes with packing frac- 
tion. This question is especially important for granular 
systems with glassy behavior, where a number of groups 
have studied the length scale related to spatially hetero- 
geneous dynamics ^ iTj . In this paper we demonstrate the 
existence of both correlations and anti-correlation in the 
free Voronoi volumes, and we measure their spatial extent 
as a function of volume fraction. 



Experiment. — Experiments are performed in a two- 
dimensional air fluidized bed, as sketched in fig.[Ha). The 
particles are a binary mixture of Teflon discs with diam- 
eters of c?s = 6 mm and d/ = 9 mm. They are confined 
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between two vertical glass plates (thickness 12 mm) sep- 
arated by a distance slightly larger than the thickness of 
the discs (3.86mm). The bed contains approximately 750 
discs of each size. 

The disc packing is tapped from below by air pulses 
flowing through a distributor of open-porous foam (Duocel 
40 PPI aluminium foam). Electrostatic charging is min- 
imized by grounding the distributor. The duration and 
strength of air pulses are controlled by two Waston Smith 
06B04604 mechanical pressure regulators and two pairs of 
Jefferson 2026 series electronic valves. Three sensors below 
the distributor are used to measure air pressure (Validyne 
DP15), humidity (Honeywell HIH-3610) and temperature 
(YSI 44033). Typical humidities and temperatures are 
3.8 ±0.3% and 24.8±0.4°C. 

The average packing fraction 0avg value is controlled by 
the type and duration of the air pulses. The use of different 
tapping modes (cf. fig.IHb)) enables us to vary 0avg from 
0.8175 to 0.8380, as shown in table [1] The most compact 
beds are obtained by first driving the bed to a new con- 
figuration using a primary pulse, and then following that 
by shorter secondary pulses to compactly the bed. The 
duration of the primary pulse and the strength of air flow 
are fixed, but the number and duration of the secondary 
pulses are adjustable. Looser packings are obtained using 
two airfiow pathways with different flow rates; when the 
primary pulse stops, the secondary pathway still provides 
some small flow, which slows down the settling discs. 

After each tap, the packing is allowed to relax for four 
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Table 1: Parameters of the air pulses. All the experiments 
start with the same primary pulse (3.2 bar measured at the 
regulator, 200 ms), either followed by several secondary pulses 
with the same pressure or accompanied by an extended pulse. 



0avg Extended pulse 
0.8175 0.5 bar X 1000 ms 
0.8183 0.45 bar X 1000 ms 
0.8209 
0.8218 
0.8231 
0.8256 
0.8277 
0.8315 
0.8337 
0.8356 
0.8366 
0.8380 



Secondary pulses 





3 X 15 ms 
5 X 15 ms 
10 X 15 ms 
3 X 30 ms 
7 X 30 ms 
15 X 30 ms 
20 X 30 ms 
15 X 30 ms + 15 X 15 ms 
30 X 30 ms + 30 x 15 ms 



seconds. Then an image of the packing is taken by a CCD 
camera with a Nikkor 50mm lens. Only the central region 
of the packing (252 x 192 mm^). five small disc diameters 
away from boundaries, is captured. The spatial resolution 
is 0.17 mm per pixel. For each experiment 8000 configu- 
rations are imaged. 

Image Processing. — To analyze the configuration 
we compute the Voronoi tessellation of the packing. Be- 
cause we consider a 2D system, volume and area are used 
inter-changeably in the following. In a first step the cen- 
ters and sizes of the discs are calculated with an accuracy 
of 0.1 pixel using a template correlation technique. 

There are two methods to identify Voronoi cells in a 
bidisperse system, radical tessellation [5] and navigation 
map [S]. Radical tessellation takes the boundaries of the 
cells as the collection of points whose tangents to neighbor- 
ing particles are equal length. The navigation map takes 
the cells as the collection of points closer to the surface of 
the corresponding particles than others in the system. In 
this work the navigation map is computed numerically on 
a grid of 1/16 pixel resolution. An example of the naviga- 
tion map is shown in fig. [2ja). 

Our packing is prepared under gravity. Sidewalls in- 
troduce slow convection rolls during an air flow pulse, and 
distributor inhomogeneities introduce gradients in the air- 
flow. It is therefore not surprising to find gradients of the 
packing fraction in the systen{l. The analysis is done for 
the central spatially most homogeneous region. The spa- 
tial variance of the local packing fraction (averaged over 
the whole 8000 taps) is calculated for circles of diameter 
ISds (see the white dashed circle in figure[2Ua)). Then the 
region with the smallest variance (in all cases smaller than 
0.0043) is chosen as the analysis region for the experiment. 
Additionally, the evolution of the global packing fraction 



and the geometrical contact number in that region are ex- 
amined to make sure that no segregation occurs during 
the course of the experiment. 

Results. — 

Individual Voronoi volumes and free volume distribu- 
tion. The distribution of individual Voronoi volumes v 
can be computed directly from the results of the naviga- 
tion map. An example is shown in fig. [2l^b). Because we 
do not want to distinguish between small and large discs, 
we follow Lechenault et al. |2] and use in the subsequent 
analysis the free Voronoi volume Vf = v ~ Vmin- The 
minimum volume Vmin is the volume that a grain would 
occupy in a hexagonal packing of identical discs. It equals 
^d^, where d is the diameter of the corresponding large 
{di) discs or small (dg) discs. While the respective mean 
free volume of small and large discs still differ {v^ = 0.135 
and v( = 0.158 in units of = irdl/A for 0avg = 0.8380), 
the success of the subsequent analysis justifies this step by 
hindsight. 

The free volume distributions of small and large discs 
and both types together are presented for two packing 
fractions in fig. ^d) and (e). For large volumes the de- 
cay is exponential but besides this feature none of the free 
volume distributions and volume distributions in fig.[2jb), 
(d) and (e) could be fitted reasonably with a gamma distri- 
bution. This result differs from 2D and 3D mono-disperse 
packings [TUHT^. 

The relative height of the two peaks in the free vol- 
ume distributions changes with 0avg- A recent study of 
the probability distributions of quadron volumes (an al- 
ternative tessellation introduced in [13]) showed that the 
position and height of these peaks can be traced to condi- 
tional probabilities of cell volumes at given coordination 
numbers [l3|. So the changes visible in figure [H^d) and (e) 
may also be related to a change of contact number. 

Two-point correlation of free Voronoi volume. To 
quantify the correlations between Voronoi volumes we use 
the two-point correlation function: 



{(V, -V^){Vj ~Vj)) 



(1) 



The spatial variance of the whole system is about 0.0435. 



where i, j correspond to two points in a distance L belong- 
ing to two Voronoi volumes, and Vi, Vj are the free volume 
of these Voronoi cells. (...) indicates averages over all the 
8000 different packings created by flow pulses. Vi and Vj 
are the mean free volumes at these points (computed sep- 
arately to remedy the effect of remaining small gradients) , 
and = {af + crj)/'2 is the corresponding variance. 

To obtain better statistics, 240 pairs of points with the 
same L are selected. For each pair Cij{L) is computed us- 
ing eq. [TJ In practice the pairs are selected in the following 
way: We construct a hexagonal "wheel" centered in the 
analysis region. The length of each edge of the wheel is 
set to be L (c.f. fig. HJe)). The six outer points and the 
center form 12 pairs that are separated by a distance L. 
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Fig. 2: (a) Voronoi cells computed using the navigation map for the packing in fig.[Hc) and labeled with different colors. The 
white dashed circle indicates the size of the region over which we perform further analysis for over 8000 individual configurations. 

(b) The probability distribution of the volume of individual Voronoi cells; the two peaks correspond to the two sizes of discs. 

(c) An illustration of the construction of pairs by a hexagonal wheel. The wheel vertices (dark dots) define 12 pairs of points. 
The free volume of the Voronoi cells to which the vertices belong is used for the two-point correlation measurement, (d) and 
(e) Free volume distributions of smaU discs, large discs, and both for two different packing fractions. The volume is normafized 
by the volume of small particles Vs — Tvdg/4:. The probabilities are normalized by the number of small, large, and both discs, 
respectively. 
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Fig. 3: Depending on (j), corr{L) exhibits both correlations 
and anti-correlations. Inset: the anisotropy measurement for 
(/"avg ~ 0.8380. The two point correlation of pairs placed hori- 
zontally (dark green) and vertically (brown), rather than con- 
structed from the hexagon. 
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Fig. 4: Three characteristic lengths: Lc, obtained from a lin- 
ear fit (dashed black line) to corr{L) for small L; Lmin, corre- 
sponding to the minimum of corr{L); and Lac, obtained from 
an exponential fit (green line) to corr{L) for L > Lmin. The 
inset shows the negative area Aneg of corr[L) as a function of 

of 4>avg- 



For each such pair, the free volume of the two Voronoi 
cells to which the two vertices belong are taken as Vi and 
Vj in cq. [T] Then the hexagon is translated to five differ- 
ent positions and rotated to four different angles for each 
position. However, the whole wheel stays inside the white 
circle in fig. [21[a). Then the 240 pairs arc averaged: 

corr(L) = ^^C.,(i). (2) 

The correlation function corr{L) is shown for two val- 
ues of (jiavg in fig. [3] For low (/)avg, corr{L) decays to 
zero and then fluctuates around it. Positive values of 
corr{L) indicate that the two free Voronoi volumes de- 
viate from the average in the same direction. For high 
(/)avg, corr{L) decreases to a negative minimum and then 
increases towards 0. Negative values of corr{L) character- 
ize anti-correlations: Voronoi cells at this distance deviate 
in opposite directions from the average free volume. 

We define three characteristic lengths to describe 
corr{L) (see fig. S]). First we do a linear fit to corr{L) 
for a range centered at half of the maximum of corr(L) 
with width ±l/2ds. The point where this fit crosses zero 
yields the length Lc- For measurements of corr{L) show- 
ing anti-correlation we define a second length, Lmim which 
is extracted from a local parabolic fit around the minimum 
of corr{L). The third length. Lac, is obtained from an 
exponential fit ranging from Lmin to the end of corr{L) 
(we allow for a small offset in this fit, but the magnitude 
of this offset is less than 0.0006 in all cases). 

To find the onset of the anti-correlations, the area Aneg 
where corr{L) < (cf. fig. |4|) is plotted as a function of 
(j)avg, as shown in inset of fig. 0] |15| . An extrapolation 
of a linear fit to Aneg to zero defines the threshold for 
anti-correlations, (f>Ac = 0.8277 ± 0.0005. For </> > (pAC, 
corr{L) can be both positive and negative; therefore, it 



can not be described by a power law [2]. 

The dependence of Lc, Lmin, and Lac on 0avg is shown 
in figure [SJ Lc and Lmin slowly decrease monotonically 
with 4>a.vg, while Lac grows approximately linearly. The 
finite extent of correlations and anti-correlations can also 
be seen in in fig. El which shows how the variance of the 
free volume of a cluster scales with the number N of parti- 
cles in the cluster. In the large N limit a linear relationship 
is recovered as predicted by the central limit theorem in 
the absence of correlations. 

Gravity breaks the isotropy in our experimental setup; 
this anisotropy is visible in the correlations plotted in the 
inset of figure [3l where pairs of horizontal and vertical 
points are averaged separately. While the qualitative fea- 
tures of corr{L) are independent of direction, the lower 
statistics of this analysis does not allow us to extract the 
corresponding characteristic lengths. 

We have also performed our analysis with the radical 
tessellation and have found that all features stay qualita- 
tively the same. Quantitatively there are slight differences: 
(jjAC = 0.828 ±0.002; Lc, Lmin and Lac are 6%, 12% and 
30% larger on average. 

Discussion. — This is the first observation of anti- 
correlation between Voronoi volumes in a granular system. 
One reason for this is that earlier studies of volume fluc- 
tuations focused on the cluster composed of neighboring 
particles [2j[3]. There the variance of average free volume 
per grain ct^ as a function of the number of grains in- 
cluded N was measured, and the scahng cr^ ^ N~°' with 
a < 1 indicated the existence of correlations. However, 
the scaling between cr^ and iV is a measure of the corre- 
lation integrated over the whole cluster size. In the large 
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Fig. 5: Pacicing fraction dependence of the characteristic 
lengths of correlation Lc and anti-correlation Lmin and Lac- 
The lengths are normalized using the average disc diameter 
davg, which corresponds to 1.23 — 1.28 times ds. 



N limit, the relation between N " and corr{L) is: 



c_ 

N 



L{N) 



r{L)LdL 



(3) 



where C is a constant proportional to the density of the 
cluster (7rL(A^)^ would be the size of the cluster). The rel- 
atively small contribution of anti-correlations to the inte- 
gral in cq.|3] would be hard to distinguish from experimen- 
tal noise (cf. fig. ^ . Therefore, a direct measurement of 
two-point correlation is necessary to find anti-correlation. 

In recent years a statistical mechanics approach for 
granular systems has been developed where the Hamilton 
function is replaced by a volume function [16]. Voronoi 
cells have been used to construct this function pl fTfllTS] . 
While the discovery of anti-correlation has consequences 
for such an approach, some properties like the distribu- 
tion of the local packing fraction have been shown to not 
depend on correlations [1^) . 

It is possible that not only the extent but even the ex- 
istence of positive correlation is an artifact of Voronoi tes- 
sellation. By definition the Voronoi tessellation assigns 
space to each particle in a certain 'equal' way. In the 
case of the navigation map, the edge of a Voronoi cell 
is the collection of points that have equal distance r to 
the surface of the neighboring particles. The increase 
in a Voronoi volume could roughly be seen as the in- 
crease in r, therefore the r of its neighboring Voronoi cell. 
Consequently the Voronoi tessellation itself gives rise to 
a positive correlation between neighboring Voronoi vol- 
umes. To verify whether this effect dominates the evolu- 
tion of Lc, we estimate the average size dvor of Voronoi 
cells in the system in the following way. In a homoge- 
neous system, ^avg could be written as the ratio of the 
average particle volume to the average Voronoi volume. 
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Fig. 6: The correlations have finite spatial extent, as shown by 
the variance of the average free volume of a compact clus- 
ter as a function of the number of grains included (open blue 
circles). In the absence of correlations the central limit theo- 
rem predicts a slope of one. While the slope for small clusters 
{N < 10) is 0.76, in the large A'^ (and therefore large L) limit a 
slope of one is recovered. Also the slope of of a hexagonal wheel 
"cluster" with side length L ~ 8ds (open triangles) equals one. 
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0.8380; data sets are shifted for clarity. 
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tem davg is in the range of 1.23 
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Given the bidisperse nature of our sys- 
1.28(is. This estima- 

davg/y/4>a.vg- RcSCaliug Lc with dvor 

shows that Lc decays slightly faster than explained by 
simple compaction. This could stem from the appearance 
of the anti-correlation beyond (pAC- It will be interesting 
to test this hypothesis with another way of tessellation, 
such as quadrons |13j . 

Above (fiAC the fluctuation of the volume of one grain 
causes more and more grains to be anti-corrclatcd, which 
is indicated by the increase of Lac a-nd the decrease of 
Lmin- Concerning the physical interpretation of these 
anti-correlations, we conjecture that (pAc might corre- 
spond to the onset of dilatancy. In 3D systems it is well 
established that loose granular material collapses under 
shear, while granular material denser than dilatancy on- 
set expands pUll^ . This expansion can be understood as 
competition between the grains for free volume. 

Our results show that correlations between the volumes 
of subunits depend on the specific system under consider- 
ation. For example, in froths there exist anti-correlations 
of the number of faces between neighboring cells . 
At the same time the volume of a foam cell is proportional 
to the number of faces These features combined 

indicate that the volume of neighboring cells should be 
anti-correlatcd. 

It would be interesting to know how Lac would develop 
at even higher density. However, we can not compactify 
the system above </> = 0.838 without partial segregation. 
In this case no anti-correlations occur within the segre- 
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gated patches. The only surviving length is Lc, which 
characterizes the size of the segregated patches. 

Conclusion. — Our measurements of the two-point 
correlation function for binary disc packings yield a short- 
ranged positive correlation over the whole range of pack- 
ing fractions 0.8175 — 0.8380. Above a volume fraction 
of 0.8277 ± 0.0005 we observe anti-correlation in the free 
Voronoi volumes. These anti-correlations reach a maxi- 
mum at a distance of about 3.5 small particle diameters. 
They then decay exponentially with distance, with an ex- 
ponent growing linearly with packing fraction. 

* * * 
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